Nonlinear localized modes in Glauber-Fock photonic lattices 
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We study a nonlinear Glauber-Fock lattice and the conditions for the excitation of localized structures. We 
investigate the particular linear properties of these lattices, including linear localized modes. We investigate 
numerically nonlinear modes centered in each site of the lattice. We found a strong disagreement of the general 
tendency between the stationary and the dynamical excitation thresholds, and we give a new definition based 
on both considerations. 
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Photonic crystals and waveguide arrays are one of the 
most fruitful areas of research in optics and nonlinear 
physics. During the last years these ideal experimen- 
tal and theoretical scenarios enabled the prediction and 
demonstration of many fundamental properties of the 
general area of discrete nonlinear systems [1,2]. Studies 
of periodic and non-periodic systems allowed to formu- 
late the necessary conditions for the excitation of non- 
linear localized structures (i.e., discrete solitons) in dif- 
ferent dimensions and geometries. For example, surface 
solitons have been an intense area of research [3] where 
stationary localized solutions exist only above a power 
threshold. 

The development of new experimental techniques has 
allowed a variety of interesting configurations. In partic- 
ular, we mention here the laser direct- writing approach 
using focused pulsed femtosecond laser radiation [4] . The 
possibilities range from ordered to disordered systems 
in ID and 2D settings. For example, main properties 
of chirped arrays have been theoretically [5] and exper- 
imentally [6] analyzed. Very recently, in the context of 
quantum optics in waveguide arrays, Glauber-Fock (GF) 
photonic lattices [7,8] have been studied. Authors found, 
for example, the appearance of a new family of quantum 
correlations which are absent in uniform arrays. 

In this letter, we study the excitation of nonlinear lo- 
calized states in a GF lattice. This system corresponds 
to a single-mode optical waveguide array whose coupling 
between sites increases as a square root. In the coupled- 
mode approximation, the evolution of the field ampli- 
tude is described by a discrete nonlinear Schrodinger-like 
equation [1, 2, 7, 8]: 

-i^T = Vn+iE n+1 + V n E n _ x + j\E n \ 2 E n , (1) 

where E n denotes the amplitude of the mode in the n-th 
waveguide for an array of N sites, z corresponds to the 
propagation coordinate along the waveguide, 7 to the 
nonlinear parameter and V n = to an increasing 

coupling coefficient (GF lattice). Without loss of gener- 
ality, we define V\ = 1 and 7 = 1 (focusing case). To 



characterize the solutions, we use the integral power (P) 
and the participation ratio (R) defined as 
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Fig. 1. Intensity distribution (\E^\ ) of linear modes. 

Linear modes of model (1) are given by E%(z) = 
\fn\2~ n l 2 R n {\ v j yf^i) exp(zA zy z), with H n being an Her- 
mite polynomial [7]. The eigenvalue spectrum {A^} does 
not form a compact band, but rather a set of well sep- 
arated discrete frequency levels. The system is invariant 
under the staggered-unstaggered transformation: 

E [-v+(N+l)/2] = ^n E lu+(N+l)/2] y„ e ^ (AT-l)/2}, 

with A [zy+(A r + i) /2 ] = -\-v+(N+i)/2]- Therefore, only 
half of the modes are shown in Fig.l for N = 59, follow- 
ing Ref. [8] . The spatial distribution of the linear modes 
is very inhomogeneous. For modes with frequencies close 
to zero (y ~ 30), light is extended across the array with 
a well defined main lobe centered at different positions of 
the lattice [see Fig.l]. For modes with frequencies closer 
to the band edges (y ~ 1 or 59), profiles are well local- 
ized. [For an homogenous DNLS lattice, all linear modes 
are spatially extended (R ~ 2N/3) without any main 
lobe or peak [1, 2]]. 
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As a consequence of the linear properties (mode pro- 
files and frequency gaps), nonlinear modes will exist in 
different frequency domains. For example, modes v = 1 
and 59 are the most localized ones {R « 9), both having 
a peak located around n ~ 55. These modes have no ex- 
tended tail and they are a result of the ramping in the 
coupling coefficient (for chirped systems, there is always 
at least one linear localized mode close to one of the bor- 
ders of the lattice [5,6]). Therefore, we expect to find a 
nonlinear solution bifurcating with zero power threshold 
close to this region. On the other hand, the center mode 
(A30 = 0, R « 18) has a main peak located at n = 1. 
This mode possesses a non-negligible tail which will for- 
bid the excitation of a zero-threshold fundamental non- 
linear solution (one-peaked bright solitons, i.e. an odd 
mode [2]). In addition, modes v = 28, 29, 31 and 32 also 
possess a main peak located at n = 1. Therefore, we ex- 
pect that different linear modes - at different frequencies 
- contribute in the formation of the profile/tail structure 
while the solution decreases its power and interact with 
the linear spectrum. For example, modes v — 26, 27, 33 
and 34 possess a main peak located at n = 2. Therefore, 
the respective nonlinear mode will not have an unique 
frequency origin, and it will probably experience some 
power threshold. By observing Fig. 1, we expect that 
the power and frequency thresholds increase as the po- 
sition of the linear mode peak increases with n (up to 
the region where we expect a zero threshold). Then, the 
power threshold should increase again because there is 
no further linear mode peaked at n > 56. 

We look for stationary solutions of Eq.(l) of the form 
En(z) = A n exp(zAz), where A n are real numbers and A 
corresponds to the solution frequency. We compute one- 
peak localized stationary solutions by using a Newton- 
Raphson iterative method starting at A ^> A59 « 14. 
Once we get a numerical solution we make a sweep in A to 
construct the whole family. For each solution, we perform 
a standard stability analysis [3] and, as a general conven- 
tion, we use black (gray) lines for stable (unstable) solu- 
tions. Fig. 2 (a) shows a Power versus Frequency diagram 
for solutions with a main peak located at n = 1. De- 
pending on the accuracy of our numerical continuation, 
we detect some "resonances" of the nonlinear solutions 
with some linear band modes (vertical dashed lines indi- 
cate different A^). If we numerically continue these solu- 
tions - by increasing the frequency - we see how the tails 
increase with an structure similar to the linear modes 
at those regions [see divergent curves in Figs.2(a)-(b)]. 
Solutions peaked close to n = 1 present similar features; 
see, for example, Fig. 2(b). Insets in Figs.2(a)-(b) show 
some profiles for these families. We see how high fre- 
quency (power) solutions are well localized and how low 
frequency solutions interact with the linear modes (see 
the tails of these profiles). If we excite solutions peaked 
at larger n, we find a phenomenology typical for surface 
modes [3]; i.e., the power and frequency decrease up to 
a point where dP/dX < 0, solutions become unstable 
and are not allowed to exist below a power threshold. 
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Fig. 2. (a) and (b) P versus propagation constant (A) for 
solutions peaked at site n = 1 and n = 5, respectively 
(vertical dashed lines correspond to different X u ). Insets 
show profiles {A n } for parameters indicated in each fig- 
ure, (c) P st versus the peak position n. 



In contrast, the solution peaked at n = 56 exhibits a 
power that drops to almost zero, because it bifurcates 
from the linear localized mode v = 59. We character- 
ize all the stationary solutions by defining the stationary 
power threshold (P s t) as the last set {A,P} before the 
solution becomes unstable [horizontal lines in Figs. 2(a)- 
(b)]. Fig. 2(c) shows P st versus the site where the solution 
has its main peak. We see an approximately square-root- 
dependence of the power increment from n = 1 up to the 
region n ~ 56, where this threshold abruptly drops to 
zero. For the further solutions, the stationary threshold 
increases again up to a maximum. 

In a next step we want to determine the dynami- 
cal power threshold for the excitation of one-site solu- 
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tions centered in different sites of the lattice. We in- 
tegrate numerically model (1) with a one-site input 
condition E n (0) = V^Ai,n > where Pq is the input 
power and no the input position. For the set {no,Po} 
we compute the space- averaged fraction of power: T = 
l/(Poz max ) f* max \E no (z)\ 2 dz, remaining at the initial 
waveguide after a fixed propagation distance z max . Fig. 3 
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Fig. 3. T as a function of Po and no- Continuous, dashed 
and dotted white lines correspond to T = 0.8, P s t and 
P ps , respectively. Inset: zoom of region n G {45,59}. 

shows our results where T = 1 (lighter color) means 
that 100% of the light is confined in the input waveguide 
whereas T = (darker color) means that all the light 
has escaped from the input position. In this figure, we 
have included a continuous white line for T = 0.8 as an 
arbitrarily chosen localization indicator. This line shows 
an increasing dynamical threshold for nonlinear states 
with an increasing position. Then, this threshold just 
decreases up to some minima at the last site of the ar- 
ray. A white dashed line indicates the stationary power 
threshold P st . The main tendency is corroborated and 
the dynamical threshold is always larger than the sta- 
tionary one. However, in the region of zero stationary 
power threshold there is an important disagreement. For 
this reason, we introduce the pseudo-stationary power 
threshold (P ps )- We define P ps as the value of the sta- 
tionary power P for which the peak amplitude is an 80% 
of the total power (| A no | 2 = 0.8P). This quantity is com- 
puted for stationary solutions and is represented with 
a white dotted line in Fig. 3. Our definition incorpo- 
rates the experimental restriction of one-site input exci- 
tations for exciting highly localized states (P ps implies 
that most of the power is contained at the peak). The 
experimental excitation of stable fundamental solutions 
like the ones sketched in insets of Figs.2(a)-(b) is not a 



trivial issue due to the long tail of the solutions. From 
Fig. 3 we see how the tendency of our definition agrees 
perfectly with the dynamical result. Nevertheless, the 
pseudo-stationary prediction will always be smaller than 
the dynamical one (for example, bulk DNLS solutions 
have a stationary threshold equal to zero and a dynam- 
ical one around 4). 

Additionally, we have studied systems with different 
functional forms of V n (exponential, quadratic, linear 
and logarithmic). All these systems possess a similar 
band structure in the sense of non-compactness and spa- 
tial distribution of linear modes. The numerical compu- 
tation of T shows similar features than the ones sketched 
in Fig. 3, except the global shape of the increment which 
is proportional to the shape of the particular V n function. 
In order to experimentally prove the threshold behavior 
and observe nonlinear localized modes in such structures, 
the range of the coupling function V n is constrained to 
the length of the sample. A large value of V n implies a 
strong coupling and a faster transport from one guide to 
another and, therefore, requires a very long crystal for 
observing such phenomenology. 

In conclusion, we studied the linear properties and the 
stationary and dynamical excitation of nonlinear local- 
ized modes in GF lattices. Contrary to homogenous sys- 
tems, we found that close to the lattice-edges solutions 
possess smaller stationary thresholds. We introduced a 
new quantity, the P ps , which showed its utility when 
comparing stationary and dynamical properties. 
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